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1. Introduction 

Nowadays, with the increasing precision of the experimental determinations of many physical 
quantities, and in some cases with the improvement of the theoretical predictions, the control over 
isospin breaking effects is becoming phenomenologically relevant. 

In the past, isospin breaking effects due to the light quarks mass difference have been ac- 
commodated within the chiral perturbation theory (ChPT) framework [1-6], while several attempts 
to compute electromagnetic effects for the hadron spectroscopy in lattice QCD have been pre- 
sented [7-10]. 

Isospin breaking effects in the Standard Model are induced both from light quark masses 
(QCD effects) and charges difference (QED effects). In the real world the separation between the 
two terms is conventional due to different additive renormalization of u, d quark masses in presence 
of electromagnetism. In the present work we consider the well-defined theoretical limit in which 
the electromagnetic interactions are switched off, and present a new method to compute the leading 
QCD isospin breaking effects with high precision. The method is based on the expansion of the 
lattice path-integral in powers of the small parameter m d — m u and is applicable in principle to any 
hadronic observable which can be computed on the lattice. 

Till now we have applied it to the computation of the kaon masses, the kaon decay constants, 
and the neutron-proton mass splitting. In the future we plan to apply the method to other physical 
quantities, to include QED corrections and to try also the calculation of next-to-leading corrections 
such as the n + -n Q mass difference. 

The results presented in this talk have been published in [1 1]. Here we present a shorter version 
of the published work: we marked in the text the parts which are discussed in more details in the 
paper. 



2. Description of the method 

Let us start by considering the evaluation of a generic euclidean correlation function {&) used 
to extract information about physical quantities as masses, decay constants, form factors etc., 

(0} = ^jD<j)0e- s , 2? = jD$e- s . (2.1) 

We can split the Lagrangian into SU (2)y symmetric and isosping violating terms, 

-Sf = Jz? - Am K(i i^ = S£un + m ud qq - Am ud qZ 3 q, (2.2) 

where q T = (u,d), m ud = (m d + m u )/2 and Am ui i = (m d — m u )/2. By expanding at first order the 
exponential of the action, S = L^-if (x), with respect to Am„j we obtain: 



J Dtp 0(l+Am ud S)e- So _ (0 Q +Am 



ud 



jD(j> (\+Am ud S)e- s o l+Am ud {S) 



o + Am ud (0S) Q , (2.3) 



where (-)o represent the vacuum expectation value in the isospin symmetric theory and S is the 
isospin breaking term S = [q^q] (x) = Y x [" M ~~ dd] (x). The correction in the denominator van- 
ishes, (5)o = 0, because of isospin symmetry. We can now describe a general recipe to be used in 
order to compute leading QCD isospin breaking effects on the lattice: 
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• consider a given correlation function in the full theory, i.e. with m u ^ rrid, and for each gauge 
configuration draw all the fermionic Wick contractions; 

• expand the up and down quark propagators with respect to Am u d according to 

G„/d(xi,X2) = GfXxi,x 2 )±Am ud J^G £ (xi,y) G £ (y,x 2 )-\ ; (2.4) 

y 

• retain the terms linear in Am Ui i and compute the corresponding diagrams. 
Eq. (2.4) can be represented diagrammatically as follows 

= > ±— <$— + ■■■, (2.5) 



where plus holds for up (light blue) and minus for down (green) quarks. Black lines from x to y refer 
to G((x — y) = (£(x)£(y)}, the propagator with the symmetric mass m u d in the isospin symmetric 
theory, whereas the cross represents the insertion of the renormalization group invariant quantity 

® = Am ud Z z £(z)£(z). 

Note that the square of the Gt propagator entering eq. (2.4) can be easily calculated on the 
lattice by using Gt itself as the source vector of a new inversion. 

3. Kaon masses and decay constants 

In this section we discuss in detail the strategy used to derive the isospin corrections to the 
kaon masses and decay constants. To this end we start by expanding the Euclidean two-point 
pseudoscalar correlation functions of kaons around the isospin symmetric point: 

s 

C K +K-(t) = Y, e ~ iP '*( a Ks(x) s Y5 u(0)) = -<^> = -<^>-<^>+^(Am Hrf ) 2 , 

v u 

s 

* d 
The spectral decomposition of C K o K o (the analysis of C K + K - proceeds along similar lines) is 

q2 

C&go&t) = Y,e- rp ' x (d~Y 5 s(x,t)s 75 d(0)} = -f- e~ E ^ + • • • , (3.2) 

By differentiating eq. (3.2) with respect to Am it is easy to see that from the ratio of the two 
correlators 
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Figure 1: Left panel: extraction of meson energies from the effective mass of Ckk{p,i). Right panel: fits of 
SCKK{P,t)/aAm^ d : as it can be seen numerical data follow theoretical expectations. The data correspond to /3 = 3.9, 
am ud = 0-0064, am^ = 0.0177, with am^ the bare mass for quark q. 



it is possible to extract the leading QCD isospin breaking corrections to kaon energies and decay 
constants. Indeed AEk appears directly in the previous equation as the "slope" with respect to t 
whereas 8Fk can be extracted from the "intercept" according to 



(m s + m ud )— j, 
M 2 K 



8F K 



Am ud 
m s + m ud 



+ 8G K -28M K . 



(3.4) 



On a lattice of finite time extent T with quark fields satisfying anti-periodic boundary conditions 
along the time direction and given our choice of the kaon source and sink operators, the pseu- 
doscalar densities, eq. (3.3) has to be properly modified. 

In this work we have used the Nf = 2 dynamical gauge ensambles generated and made publicly 
available by the European Twisted Mass Collaboration. These gauge configurations have been 
generated by using the so called Twisted Mass lattice discretization of the QCD action [12]. For the 
present analysis we have used the same 13 gauge ensembles taken into account in [13], using the 
same values of renormalization constants and lattice spacings there mentioned, and considering a 
statistics of about 150 configurations per ensemble. The computation of all the correlation functions 
has been carried out on AURORA machine in Trento. The full set of parameters can be found in 
the appendix A of [11], where we discuss also the relevance of isospin breaking effects at finite 
lattice spacing caused by the Twisted Mass regularization. 

As can be seen from Figure 1, 8Ckk(p,i) is determined with high precision, given the strong 
statistical correlation existing between the numerator and the denominator of the ratio in eq. (3.3). 

Our results do not show a visible dependence with respect to m ud within the quoted errors. 
Therefore we simply extrapolated the correction to the meson mass square (which is a finite quan- 
tity in the chiral limit) to the continuum using the expression: 



AMf 
Am ud 



{m ud ,a) 



AMI 
Amud 



QCD 



+ Cmct 



(3.5) 



In the left panel of Figure 2 we show the continuum extrapolation of AM\/ Am ud . Adding a linear 
and a (chiral) logarithmic term according to the formulae obtained within the unitary SU(3)l x 
SU (3)r in ref. [1] leads to compatible results for [AMg/Am ud ]^ CD , with a 3% difference which we 
considered as an estimate of systematic error, adding it in quadrature to the lattice uncertainty. 
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Figure 2: Left panel: combined chiral and continuum extrapolations of AM^/Am ud , in physical units. Right panel: 
combined chiral and continuum extrapolations of Si^/Am^. Black points correspond to the coarser lattice spacing, 
a — 0.098 fm, dark magenta points correspond to a = 0.085 fm, green points to a = 0.067 fm and blue points to 
a = 0.054 fm. Red lines are the results of the continuum extrapolations. 

In the case of dF^/Amud (right panel of Figure 2), the dependence upon m ud is significant, and 
we have included in the fitting function the leading and next-to-leading terms expanded in powers 
of m lid /m s plus a lattice artifact term, i.e. 

qcd / \ 

+ C F a 2 + Si [m ud - m^ D \ + B 2 m ud log 



' 8F K - 




' 8F K - 


Am ltd _ 


(mud, a) = 


Am ud 



mud 



QCD 

ud , 



(3.6) 



The systematics associated to this extrapolation has been estimated by replacing the logarithmic 
term with a quadratic one, and it has been found of the order of 5%. 

Using the value m^ D = m^ D (MS,2GeV) = 3.6(2) MeV from refs. [13, 14] we obtain 

QCD 



AM Z K 
Ain ud 



■ 2.57(8) x 10 j MeV, 



MS,2GeV 



8F K 
Antud 



QCD 



MS,2GeV 



3.3(3) x HT 3 MeV -1 . (3.7) 



-l 



Having neglected the QED effects from our calculations, we cannot directly use the exper- 
imental determination of M^ Q —M^+ to extract [m d — m„]@ CD . By taking the Chiral Perturbation 
Theory (ChPT) estimate of the electromagnetic corrections to M^ —Mg + we can obtain the theoret- 
ical value of the kaons squared mass difference in absence of electromagnetic interactions [7-10], 
ref. [14], which reads: 

]QCD 



[M 2 Ka -M 2 K+ ]^ u = 6.05(63) x 10 3 MeV 2 . 



(3.8) 



Using eq. (3.8) with our numerical results given in eq. (3.7), we get 



m d - m u ] QCD (MS,2GeV) = 2.35(8)(24) MeV 



F K o - F K 
F K 



QCD 



0.0078(7) (4) (3.9) 



where the first error comes from our calculation and combines in quadrature statistics and system- 
atics while the second comes from the uncertainty on QED effects. At first order in Am ud , due to 
the fact that pions don't get corrections and that K + and K° get opposite corrections, we have 

QCD 

= -0.0039(3) (2). (3.10) 



Fk/Fji 

This is significatively higher than the estimate obtained in ref. [6] by using chiral perturbation 
theory, where the value -0.0022(6) was found. 
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Figure 3: Left panel: Correlation functions 5Cmv(0 I a ^ m ud- R'ght panel: Chiral extrapolation of AM^/Am ll( {. The 
data are at fixed lattice spacing a = 0.085 fm for different values of m ue j. 

4. Nucleon masses 

Having determined Am^ D , we can now predict the QCD part of the difference between the 
masses of the neutron and of the proton. We consider for proton the correlation function 



e abc {u a Cy 5 dl)u, 



l±y° 



1 ±y° 

£def — ^ — u d [u e Cy 5 d f ) 



(0)>. 



and similarly for neutron, with «fj>(/. To increase statistic, we consider the combinations 

C nn (t) = C+ (t) - C- n (T - 1) , C pp (t) = C+Jt) - C-JT - 1) . 



(4.1) 



(4.2) 



Quark fields have been "Gaussian smeared" according to the values found in ref. [16]. The cor- 
rection to such correlation function with respect to the isospin symmetric theory can be found 
by expanding eq. (4.2) and computing all the required contractions, as already shown in previous 
section in the case of kaons. 

The extraction of physical information from nucleon euclidean two point functions proceeds 
along the same lines described in detail in the case of the kaons. By extracting the slope in 
t of 5Caw(?), we can determine AM^ = (M n —M p )/2. In the left panel of Figure 3 we show 
SCw/v^/oAm^, that we have fitted with the linear function: 8Cnn = c + tAMpj ■ 

In the right panel of Figure 3 we show the chiral extrapolation of AMm/Am^ performed by 
using the following fitting function [17] 



AM N 
Am ud 



{mud) 



AM N 
Am lld 



QCD 



+ B N (m ud -m^ D 



(4.3) 



By using the results of the fit and the value of Ampf given in eq. (3.9), we get 



[M n -M p 



QCD = 2Am% D 



AMa, 



Am 



ml 



QCD 



2.8(6)(3)MeV 



\QCD 



6.05 x 10 3 MeV 



2 ' 



(4.4) 



where the first error takes into account the lattice uncertainties while the second comes from the 
uncertainty on QED contribution. This is our best estimate at present but it has been obtained at 
fixed lattice spacing and with limited statistics. A refinement of this calculation is in progress. 
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5. Conclusions and Outlooks 

In this talk we have proposed a new method to compute with high precision the QCD isospin 
breaking effects in relevant physical quantities at the lowest non trivial order in the up-down mass 
difference. The method can be easily extended with minor modifications to higher orders. We 
have computed the corrections to meson and nucleon masses and meson decay constants, showing 
that, in spite of the limited statistics, our approach is already competitive, or even better, than other 
calculations based on the effective QCD chiral lagrangian. 

To obtain the complete physical results, our method has to be combined with the calculations 
of the electromagnetic corrections, which is currently under investigation. In this work, for a com- 
parison with calculations in different theoretical frameworks, we have taken the electromagnetic 
corrections to the meson masses evaluated in ref. [14]. 

As the method looks very promising, we are planning to extend this work to other physical 
observables, such as the form factors of semileptonic decays, for which we have presented 
preliminary results in the full paper [1 1]. 
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